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ABSTRACT 

This  study  seeks  to  address  problems  in  manufacturing,  and  in  particular,  schedule  planning.  Black  Widow  Optimization 
(BWO)  metaheuristic  is  the  algorithm  employed  ,  which  can  identify  the  optimal  usage  of  the  tariables  in  the 
manufacturing  process.  Computational  experiments  are  carried  out  with  data  prorided  by  an  unspecijred  manufacturer, 
under  varying  constraints.  The  jindings  indicate  that  various  problems  can  be  solved  by  BWO  with  superior  efficiency  to 
alternative  approaches  and  within  a  reasonable  timeframe.  The  findings  also  suggest  that  future  studies  should  seek  to 
provide  a  comprehensive  comparison  with  alternative  metaheuristics,  and  more  complex  problems  should  be  addressed, 
such  as  those  involving  different  product  groups  or  demand patterns. 
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1.  INTRODUCTION 

There  are  two  basic  types  of  optimization  algorithm;  approximation  optimization  algorithms  and  conventional 
optimization  algorithms.  The  former  can  be  employed  to  find  the  best  value  for  a  solution  when  a  conventional 
optimization  algorithm  is  used.  However,  while  it  is  effective  with  small  problems,  it  can  be  much  harder  to  find  a 
solution  in  the  case  of  large  problems.  Notably,  where  problems  are  of  the  NP-Hard  type,  a  slight  increase  in  the 
size  of  a  problem  results  in  a  much  longer  time,  which  is  necessary  to  determine  a  solution.  In  this  scenario,  the 
approximation  principle  offers  the  best  approach  to  determine  a  direct  solution.  The  techniques  of  this  type  can  be 
used  to  solve  much  larger  and  more  complex  problems.  The  branch  of  metaheuristics  deals  with  using  the 
approximation  principle  to  discover  the  optimal  solution.  It  allows  complex  problems  to  be  processed  rapidly,  with 
the  traveling  salesman  problem  among  the  common  examples.  The  different  metaheuristic  approaches  can  be 
separated  into  four  categories.  The  first  type  is  known  as  evolutionary  algorithms,  such  as  Evolutionary  Strategy 
(ES)  (Cheng  et  al.,  2007),  and  Genetic  Algorithm  (GA)  (Cus  and  Balic,  2003).  Second  group  is  Animal  Algorithm. 
The  popular  member  in  this  group  are  Grey  Wolf  Optimization  algorithm  (GWO)  (Mirjalili  et  al.,  2014),  and 
Artificial  Bee  Colony  (ABC)  (Karaboga,  2005).  Physical  algorithm  is  a  third  group.  This  group  includes  Simulated 
Annealing  (SA)  (Granville  et  al.,  1994),  Gravitational  Search  Algorithm  (GSA)  (Rashedi  et  al.,  2009),  and  Thermal 
Exchange  Optimization  (TEO)  (Kaveh  and  Dadras,  2017).  The  last  group  simulates  Human  Behavior  such  as  Tabu 
Search  (TS)  (Lokketangen  et  al.,  1994)  and  Teaching-Learning-Base  Optimization  (TLBO)  (Rao  et  al.,  2011). 
These  techniques  generate  solutions  which  might  not  represent  the  optimal  answer,  but  instead  produce  a  group  of 
answers.  The  working  process  continues  to  operate  until  the  conditions,  which  were  specified  from  the  outset  are 
achieved.  Typical  examples  of  this  approach  include  the  TEO,  and  GA  which  are  techniques  that  involves 
adjustments  whereby  different  methods  can  be  combined  in  order  to  eliminate  any  shortcomings  and  find  new  ways 
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to  reach  a  solution  through  a  combination  of  different  principles.  Each  technique  will  find  answers  in  its  own  specific  way. 
This  research  paper  describes  an  approach  to  discovering  optimal  solutions.  This  is  known  as  BWO  (Hayyolalam  and 
Kazem,  2020),  which  is  a  type  of  metaheuristic  algorithm  initially  introduced  in  2020  by  Hayyolalam  and  Kazem. 

2.  BLACK  WIDOW  OPTIMIZATION  ALGORITHM  (BWO) 


Figure  1:  FIowchart  of  the  BWO  (Hayyolalam  and  Kazem,  2020). 

The  riowchart  for  the  algorithm  proposed  by  this  study  appears  in  Figure  1.  The  algorithm  begins  with  a 
population  of  spiders  whereby  each  of  the  spiders  can  be  understood  to  represent  one  possible  answer.  The  spiders  then 
pair  up  to  reproduce,  creating  a  new  generation.  Female  black  widows  are  known  to  devour  the  male  during  mating  or 
upon  the  conclusion  of  mating,  as  can  be  seen  in  Figure  2.  The  female  then  stores  the  sperm  in  her  cae  before  they  are 
allowed  to  enter  the  egg  sac,  as  can  be  observed  in  Figure  3(a).  The  baby  spiders  appear  after  11  days,  and  will  live 
together  on  the  web  of  the  mother  spider  for  up  to  a  week,  often  eating  each  other,  as  shown  in  Figure  3(b).  They 
eventually  leave  the  web  when  they  are  carried  away  by  the  wind. 


Figure  2:  A  Female  Black  Widow  Spider  in  Her  Web 
(Hayyolalam  and  Kazem,  2020). 
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2.1  Initial  Population 


In  the  context  of  an  /V.a|  -dimensional  problem  of  optimization  problem,  one  widow  represents  an  array  of  1  x  /Vvar 
which  can  be  seen  as  a  problem  solution.  The  array  has  the  definition  given  as: 


Widow=  \XVX2,...,TLN  ] 

1  ^  1  y  var 


(1) 


Each  variable  (x, ,  X7 , X  ,,  )  serves  as  a  tloating-point  number.  The  widow  fitness  can  be  determined  by 

1  'Vvar 

assessina  the  fitness  function  /  ateach  widowgiven  as  (x, ,  X, , X  v.  )  .  Therefore, 

1  z  yvvar 

Fitness  =  /  (widow)  =  f(xx,X2,...,  XN^  )  ,  (2) 

In  order  for  the  optimization  algorithm  to  begin,  it  is  necessary  for  one  candidate  widow  matrix  of  the  size 
N  x  /Vva|.  to  be  created  with  its  starting  spider  population. 

2.2  Procreate 

In  order  to  achieve  reproduction,  the  algorithm  makes  use  of  an  array  known  as  alpha,  which  is  developed  as  the  widow 
array  using  random  numbers,  where  the  resulting  offspring  are  generated  through  the  use  of  Ot  with  the  equation  given 

below  (Equation  3)  where  Xx  and  Xn  represent  the  parents  while  yx  and  yn  are  the  offspring. 

fy,  =  axxj  +(l-a)xx2 

[y2  =axx2+(l-a)xXl 

This  process  is  carried  out  repeatedly  for  a  total  of  /Vvar/ 2  times,  although  there  should  be  no  duplication  of 
numbers  selected  at  random.  At  the  end,  the  offspring  and  mother  can  be  added  to  an  array  before  undergoing  sorting  on 


www.tivrc.ore 


SCOPUS  Indexed  Journal 


editor@tjprc.  org 


1412 


Somboon  Sukpancharoen,  Surachat  Chantarachit  &  Udomsak  Jantontapo 


the  basis  of  their  fitness  values. 

2.3  Cannibalism 

There  are  three  different  types  of  cannibalism  that  occur.  Sexual  cannibalism  arises  when  the  female  eats  the  male  either 
during  or  after  mating.  The  algorithm  is  able  to  distinguish  between  male  and  female  spiders  by  considering  their  fitness 
values.  Sibling  cannibalism  occurs  when  stronger  baby  spiders  eat  weaker  baby  spiders.  The  algorithm  applies  a 
Cannibalism  Rating,  or  CR,  which  sets  the  number  of  spiders  which  will  survive.  It  is  also  possible  for  baby  spiders  to  eat 
the  mother  spider.  The  algorithm  uses  the  fitness  value  to  establish  whether  baby  spiders  are  weak  or  strong. 

2.4  Mutation 

At  this  point,  the  random  selection  of  Mutepop  number  of  individuals  from  the  population  takes  place.  The 
solutions  chosen  will  randomly  have  two  of  their  elements  from  the  array  rearranged,  as  Figure  4  shows.  In  order  to 
calculate  Mutepop,  the  mutation  rate  is  used. 


X, 

x: 

Y 

/xnvar 

2.5  Convergence 


As  would  be  the  case  for  other  evolutionary  algorithms,  there  are  three  stop  conditions  which  can  be  applied:  (a)  a 
predefined  number  of  iterations.  (B)  Observance  of  no  change  in  the  fitness  value  of  the  best  widow  for  several  iterations. 
(C)  Reaching  to  the  specified  level  of  accuracy. 


2.6  Parameter  Setting 

In  the  BWO  algorithm  under  examination,  certain  parameters  are  required  in  order  to  achieve  improved  results.  These 
particular  parameters  are  the  procreating  rate  (PP),  cannibalism  rate  (CR),  and  mutation  rate  (PM).  The  values  selected  in 
this  research  for  these  particular  parameters  can  be  seen  in  Table  1. 


Table  1:  The  Set  Parameter  Yalues 


Parameter 

Value 

PP  =  procreate  rate 

0.70 

CR  =  cannibalism  rate 

0.45 

PM  =  mutation  rate 

0.60 

3.  MANUF ACTURIN G  MODEL 

3.1  Manufacturing  Model  1 

Model  1  allows  the  analysis  of  an  experimental  design  in  order  to  determine  the  influence  of  the  cutting  parameters 
including  feed  rate,  cutting  speed  and  depth  of  cut  on  surface  roughness  in  turning  process  (Pansare  and  Kavade,  2012). 
The  model  below  was  proposed  by: 

Ra  =8.11-0.0217A-25.9£-6.37C+0.0563AB+0.0153AD  +  19.4£C  (4) 

where  A  is  the  cutting  speed  (m/min),  B  is  the  feed  rate  (mm/rev),  C  is  the  depth  of  cut  (mm),  and  Ra  is  the 
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surface  roughness  (|jm).  The  formulation  of  the  turning  optimization  problem  takes  place  as  shown  below: 

Minvmze  Ra  =f(A,B,C) 

Subject  to: 

150<A<250  (m/min);0.1<5<0.2  (mm/rev);0.5<C<1.5  (mm)  (5) 

3.2  Manufacturing  Model  2 

The  objective  of  the  next  model  is  to  discover  the  optimal  values  for  the  parameters  in  turning  operations  to  maximize  the 
MRR  (Material  Removal  Rate)  (Shivakoti  et  al.,  2012).  To  calculate  MRR,  the  following  regression  equation  was 
developed: 


MRR  =  1 .42  - 1 .83A  -  0.9S  +  10C  + 103  AB  - 1 12AC  +  0.0000 14 BC 


(6) 


where  A  is  the  feed  rate  (mm/r ev),  B  is  the  spindle  speed  (rpm),  and  C  is  the  cutting  speed  (m/min).  Equation  7 
was  chosen  as  the  objective  function  via  formulation  of  the  optimization  problem  to  determine  the  most  suitable  cutting 
parameters  as  follows: 

Maximize  MRR  =  f(A,B,C) 

Subject  to: 

0.62  <A<  0.98  (mm/rev);  40  <  5  <  1000  (rpm);  3.5  <C  <95.5  (m/min)  (7) 

3.3  Manufacturing  Model  3 

The  goal  is  to  apply  the  regression  equations  in  order  to  generate  mathematical  models  (Saravanakumar  et  al.,  2012)  of  the 
Material  Removal  Rate  (MRR)  and  surface  roughness  Ra  in  the  following: 


MRR  =  19158-  298A  - 1 121 36B  +  9 1493C  +  1749AS  +  1417AC  +  537343SC  -  7880ASC 
R  =  23 .6  -  0.33 1A  - 1 10S  -  88C  + 1 66  AS  + 1 .29AC  +  463SC  -  6.93  ASC 


(8) 


(9) 


where  A  is  the  cutting  speed  (m/min),  B  is  the  feed  rate  (mm/rev),  and  C  is  the  depth  of  cut  (mm).  To  identify 
optimal  cutting  parameters,  Equation  8  must  be  maximized  while  Equation  9  must  be  minimized  with  regard  to  the 
machining  parameter  constraints  given  as: 


60  <  A  <  80  (m/min);  0. 15  <  S  <  0.25  (mm/rev);  0. 1  <  C  <  0.25  (mm) 

3.4  Manufacturing  Model  4 


(10) 


In  this  problem,  there  were  five  machines  and  20  products.  Table  2  shows  the  time  taken  for  the  processing  of  each 
product.  The  precise  solutions  or  product  sequences  in  this  problem  are  shown  as  follows: 

Table  2:  Product  Processing  Time  for  Each  Machine  for  the  Manufacturing  Model 


Product 

Processing 

Time  C  (mins)Machine 

Due  Time 
d  (mins) 

Penalty  pt 
($/min) 

1 

2 

3 

4 

5 

A 

62 

77 

65 

79 

81 

850 

55.00 
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B 

47 

44 

84 

50 

67 

1,210 

35.00 

C 

97 

77 

99 

33 

26 

350 

15.00 

D 

87 

2 

88 

7 

91 

575 

90.00 

E 

38 

62 

95 

43 

15 

810 

100.00 

F 

67 

10 

86 

11 

39 

1,060 

90.00 

G 

31 

8 

69 

61 

7 

475 

30.00 

H 

24 

64 

30 

41 

25 

655 

30.00 

I 

72 

37 

84 

66 

14 

555 

25.00 

J 

37 

77 

60 

75 

82 

uoo 

80.00 

K 

47 

80 

91 

6 

17 

890 

80.00 

L 

49 

86 

6 

89 

92 

1,055 

45.00 

M 

25 

5 

8 

42 

6 

815 

85.00 

N 

48 

8 

55 

82 

7 

700 

40.00 

O 

97 

76 

54 

64 

66 

1,205 

40.00 

P 

16 

27 

11 

5 

64 

910 

70.00 

Q 

43 

69 

61 

6 

85 

360 

75.00 

R 

41 

5 

28 

32 

38 

715 

90.00 

S 

44 

26 

7 

19 

40 

480 

60.00 

T 

10 

37 

23 

29 

100 

850 

15.00 

4.  RESULTS  AND  DISCUSSIONS 

For  the  purposes  of  this  study,  BWO  was  written  using  a  MATLAB  2019,  a  computer  program  and  the  computational 
experiments  were  conducted  using  HP  laptop.  Assessments  were  then  made  of  the  BWO,  TEO,  and  GA  so  as  to  compare 
the  overall  performance.  For  each  algorithm,  the  main  parameters  have  their  own  intluence  over  performance  in  achieving 
the  optimal  quality  of  solution  or  shortest  computational  time.  The  best  parameter  choices  for  the  problems  are  determined 
by  considering  previous  experiments,  with  each  algorithm  using  values  suggested  in  the  literature.  Development  of  the 
parameter  values  was  accomplished  through  the  experimental  designs,  whereupon  the  hndings  were  investigated  to 
determine  the  quality  of  the  solutions.  For  the  final  BWO  algorithm,  the  parameters  were  set  in  the  form,  as  shown  in  Table 
2.  Each  of  the  algorithms  then  underwent  execution  using  the  maximum  evaluation  number  up  to  10,000,  and  for  each 
problem,  the  algorithm  performed  30  repetitions. 

The  performance  of  a  number  of  different  metaheuristics  is  listed  in  Tables  3,  4,  and  5,  in  the  context  of  roughness 
minimization,  maximization  of  the  material  removal  rate,  and  optimization  problems  with  various  constraints  and 
limitations.  The  study  takes  into  consideration  four  separate  measures  of  performance  for  comparison  purposes.  These 
include  the  average,  maximum,  minimum,  and  standard  deviation  for  both  the  execution  time  and  the  yield.  The  BWO  was 
shown  to  be  capable  of  approaching  the  optimal  point  with  superior  performance  to  the  alternatives.  One  of  the  key 
strengths  of  BWO  is  its  efficiency  of  computation,  as  shown  in  Figure  6.  In  the  case  of  an  earlier  optimization  setting,  all  of 
the  tested  metaheuristics  were  capable  of  delivering  results  close  to  the  optimum.  For  most  of  the  techniques,  the  standard 
deviation  is  very  small,  but  in  the  case  of  BWO,  it  is  very  near  to  zero.  In  the  case  of  GA,  the  time  taken  for  computation  is 
lower  than  for  all  the  alternatives. 


Table  3:  Experimental  Result  on  the  Manufacturing  Model  1 


Model  1 

BWO 

TEO 

GA 

Yield 

Time 

Yield 

Time 

Yield 

Time 

Average 

0.036 

69.9 

0.085 

70.6 

0.06 

69.8 

Std.  Dev. 

0.004 

1.25 

0.05 

1.13 

0.08 

1.01 

Max 

0.038 

71.87 

0.102 

71.89 

0.12 

71.84 

Min 

0.032 

68.42 

0.039 

68.38 

0.04 

68.56 
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Table  4:  Experimental  Result  on  the  Manufacturing  Model  2 


Model  2 

BWO 

TEO 

GA 

Yield 

Time 

Yield 

Time 

Yield 

Time 

Average 

99,668 

65.5 

99,651 

66.0 

99,656 

60.7 

Std.  Dev. 

23.7 

1.0 

41.43 

1.0 

40.02 

0.9 

Max 

99,686 

67.5 

99,685 

67.6 

99,685 

63.4 

Min 

99,651 

64.0 

99,628 

64.3 

99,628 

58.2 

Table  5:  Experimental  Result  on  the  Manufacturing  Model  3 


Model  3 

BWO 

TEO 

GA 

Yield 

Time 

Yield 

Time 

Yield 

Time 

Average 

2,003 

66.0 

2,028 

66.1 

2,024 

60.7 

Std.  Dev. 

88.5 

1.04 

72.7 

1.15 

96.7 

1.02 

Max 

2,066 

67.5 

2,080 

68.0 

2,093 

63.2 

Min 

1,940 

64.4 

1,977 

64.3 

1,956 

58.2 

This  was  a  bigger  and  more  complex  problem.  In  considering  the  details  of  the  five  machines  and  20  products  in  a 
given  factory,  it  is  apparent  that  BWO  offers  the  smallest  penalty  at  US  $104,020,  while  the  order  of  production  is  [L  J  Q  E 
FNPDARCMHTKGIBOS].  The  next  and  ftnal  approaches  are  GA  and  TEO,  which  offer  the  next  smallest 
penalty  values  at  US  $107,480  for  GA  and  US  $1 10,365  for  TEO.  These  results  can  be  observed  in  Table  6. 


Table  6:  The  Findings  for  the  Minimum  Make  Span  in  the  Scheduling  Problem 


Methods 

Sequences 

Minimum  make  span  (US$) 

BWO 

[LJQEFNPDARCMHTKGIBOS] 

104,020 

TEO 

[JQLEFPNCTDARHMKGIBOS] 

110,365 

GA 

[JQLETFPNDARCHMKGIBOS] 

107,480 

The  results  of  the  statistical  analysis  of  the  performance  of  BWO  and  the  alternative  algorithms,  GA  and  TEO,  can  be  seen 
in  Table  7.  The  table  shows  that  BWO  has  a  better  mean  value  than  the  alternatives,  but  the  robustness  of  TEO  was  shown 
to  be  superior  to  BWO,  as  indicated  by  the  standard  deviation  (Std.  Dev.).  The  CPU  times  which  are  given  in  seconds  do 
not  show  any  major  differences  between  the  three  approaches. 


Table  7:  A  Comparison  of  the  Data  Analysis  Findings  for  the  Different  Algorithms 
Addressing  the  Scheduling  Problem 


Methods 

Mean 

Best 

Worst 

Std.  Dev. 

CPU.  time  (s) 

BWO 

104,320 

104,020 

104,595 

420.2 

74.9 

TEO 

110,314 

110,365 

110,410 

28.4 

87.6 

GA 

107,951 

107,480 

108,435 

675.2 

76.0 

It  is  necessary  to  start  with  a  fixed  fitness  constraint,  if  the  three  methods  are  to  be  applied  and  tested  fairly  in 
order  to  present  the  correlation  between  the  objective  function  value  and  the  number  of  iterations.  Figure  5  presents  the  rate 
of  convergence  at  population  size  set  of  50.  The  data  show  that  GA  offers  much  faster  rates  of  convergence  than  the 
alternative  approaches,  with  the  drawback  that  the  GA  solution  is  unable  to  move  beyond  the  local  optimum  when 
attempting  to  solve  the  scheduling  problem.  The  BWO  algorithm  provides  a  slower  convergence  rate  than  the  GA 
algorithm,  but  does  reach  a  global  optimum.  Furthermore,  the  answer  generated  after  testing  with  the  number  of 
evaluations  set  at  10,000  revealed  that  the  BWO  approach  was  superior  to  the  two  alternatives. 
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Figure  5:  Comparison  of  Convergence  Rates  for  the  Pressure  Scheduling  Problem. 

5.  CONCLUSIONS 

This  research  study  assesses  the  effectiveness  of  the  novel  BWO  metaheuristic  algorithm  through  testing  in  the  context  of 
different  manufacturing  problems  where  both  maximization  and  minimization  methods  are  used,  in  addition  to  solving 
scheduling  problems.  The  data  reveal  that  BWO  is  able  to  compete  effectively  with  other  approaches,  and  delivers  a  higher 
degree  of  accuracy  than  the  alternatives.  The  GA  produced  the  best  performance  in  the  context  of  computational  effort, 
although  the  yields  were  not  as  good  as  others.  The  findings  in  this  study  suggest  that  future  work  should  involve  a  more 
complete  comparison  with  the  alternative  metaheuristics,  and  more  complex  problems  should  be  addressed,  such  as 
product  groups  and  differing  demand  patterns.  Moreover,  it  is  known  that  further  research  into  parameter  values  usually 
lead  to  enhanced  solutions  when  the  final  generation  to  the  optimal  solution  vector  is  reached. 
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